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ABSTRACT 

We present an analysis of high-resolution synoptic spectra of ten F- and G-type stars, 
seven of which exhibit periodic radial velocity variations due to the presence of one or 
more substellar companions. We searched for subtle periodic variations in photospheric 
line asymmetry, as characterized by line bisectors. In principle, periodic variations in 
line asymmetry observed at lower spectral resolution could mimic the radial velocity 
signature of a companion, but we find no significant evidence of such behavior in our 
data. Observations were obtained from 1998 March to 1999 February using the National 
Solar Observatory (NSO) 1.52-m McMath-Pierce Solar Telescope Facility on Kitt Peak 
in conjunction with the solar-stellar spectrograph, achieving a resolving power of 1.2 X 
10 5 . To characterize line asymmetry, we first measured line bisectors for the unblended 
Fe I photospheric line at 625.26 nm. To improve sensitivity to small fluctuations, we 
then combined points in each bisector to form a velocity displacement with respect 
to the line core. We searched for periodic variations in this displacement, finding no 
substantial difference between stars with substellar companions and those without re- 
ported companions. We find no correlation between bisector velocity displacement and 
the known orbital phase of substellar companions around our target stars. Simulations 
of a periodic signal with noise levels that mimic our measurement errors suggest that 
we can exclude bisector variations with amplitudes greater than about 20 m s _1 . These 
results support the conclusion that extrasolar planets best explain the observed periodic 
variations in radial velocity. 
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1. Introduction 

High precision radial velocity studies of nearby stars have resulted in the first detections of 
extrasolar planetary systems, consisting primarily of “massive Jupiters” (Mayor & Queloz 1995; 
Marcy & Butler 1998 and references therein). These detections rely on achieving precisions of 
3-15 ms -1 . Stellar activity in the form of spots and plages, non-radial pulsations, and convective 
inhomogeneities combined with rotation can distort stellar photospheric lines beyond the intrinsic 
asymmetries arising from surface convection, thereby contributing error in the form of “stellar 
noise” to sensitive radial velocity measurements (Saar & Donahue 1997; Hatzes 1996). Further 
improvement in radial velocity precision beyond this intrinsic variability limit will require increased 
attention to line profile perturbations caused by stellar atmospheric phenomena. In particular, 
changes in line shape arising from stellar atmospheric motion can mimic small radial velocity 
variations at the spectral resolution ( R ~ 5-7 X 10 4 ) typically utilized for planet searches. 

Convection in the photosphere consists of rising, hotter material (associated with bright gran- 
ules on the Sun) and falling, cooler material (e.g., the intergranular lanes), yielding a superposition 
of two profiles shifted relative to each other. The net profile is asymmetric because of unequal 
contributions from the rising and falling components. In the Sun, the asymmetry appears as a 
C-shaped curvature of the line bisector, which is the locus of points midway between the red and 
blue wings of a spectral line. The bisector is defined and discussed extensively by Gray (1988 and 
references therein), especially within the context of convection in late-type stars. 

Gray (1997) originally suggested that subtle temporal variations in (R ~ 10 5 ) line bisectors 
for 51 Peg might be the source of radial velocity variations attributed to the gravitational influence 
of a planetary companion. Gray’s (1997) challenge of the planetary companion hypothesis was 
based on a sparse data-set that appeared to show regular variations in the bisector span of an 
Fe I photospheric line with a period equal to the orbital period of the putative planet. However, 
subsequent investigations by Gray (1998) and Hatzes, Cochran & Bakker (1998) based on an 
intensive series of observations of 51 Peg did not detect any systematic changes in the line shapes 
that could be misinterpreted as radial velocity shifts. Thus, the data are most compatible with the 
presence of a sub-stellar mass companion to 51 Peg. 

Following the challenge by Gray (1997), we initiated a long-term program of high spectral 
resolution observations of 7 late-type stars with reported planetary companions along with 3 addi- 
tional late-type stars with no known substellar companions down to the sensitivity limits of current 
radial velocity monitoring programs. We note that Franqois et al. (1999) also conducted an in- 
tensive series of high spectral resolution observations of 51 Peg and v And over 4 nights at the 
Canada-France-Hawaii Telescope (CFHT). These investigators did not observe any regular change 
in line bisector shape that could give rise to spurious Doppler shifts occurring with the period of 
the planetary companion. Our investigation extends these studies through the addition of a longer 
temporal baseline of relatively frequent observations of several late-type stars reported to have 
planetary companions. 
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2. Observations and Reduction 

The spectra were obtained with the 1.52-m Main telescope of the NSO/Kitt Peak McMath- 
Pierce Solar Telescope Facility in conjunction with the stellar spectrograph. The stellar spectro- 
graph is described by Smith & Giampapa (1987). We utilized the echelle grating with the 180-mm 
Nikon transfer lens, TI 800 X 800 UV-enhanced CCD and a 10 X 200-mm Bowen-Walraven image 
slicer. This configuration yielded a resolution of R = 1.20 X 10 5 as measured directly from the full 
widths at half-maximum (FWHMs) of the ThAr comparison source lines. All spectra were centered 
on the Fe I photospheric absorption line at 625.3 nm used by Gray (1998) in his study of 51 Peg. 
The program stars are presented in Table 1. 

A quartz lamp and a ThAr comparison source were utilized for flatfielding and wavelength 
calibration, respectively. A barycentric correction was applied so that all spectra for a given 
program star are on a uniform wavelength scale. While not essential for the analysis of profile 
shapes, this facilitated the combining of multiple frames, including the construction of a very high 
signal- to- noise ratio average spectrum that is used as a template spectrum for the error analysis 
(§ 2 . 2 ). 

Our initial inspection of the reduced spectra revealed line depth variations in photospheric 
features on a night-to-night basis. This apparent variability of the line strength is an artifact 
which likely is related to the uneven illumination of the slicer assembly by the stellar image in 
conditions of variable seeing. Extraction of the individual slicer spectra showed that the more 
faintly illuminated slices yielded shallower absorption features. To each spectrum, we added an 
empirical “background” correction, x, and then renormalized the continuum, obtaining 


Scorr — 


S + X 


(1) 


1 + X ’ 

where 5 is the original reduced spectrum and S corr is the corrected spectrum. This approach 
yielded individual slicer spectra with equal line depths. For 55 Cnc, the background correction 
required to match all the slices also yielded line depths in agreement with a spectrum obtained 
independently by G. Marcy using the Lick Hamilton echelle. The excellent agreement between the 
spectra validated our empirical correction scheme. 


We constructed a reference spectrum for each star by reducing and summing spectra from 
three separate images. The remaining spectra for that star were then mapped to the reference 
spectrum, using a nonlinear least-squares fitting routine that quadratically mapped each individual 
wavelength scale onto the reference wavelength scale, fit a polynomial to match the continua, and 
applied a background correction of the form in equation 1 to yield a set of consistent spectra. 
Finally, we created a high signal- to-noise (S/N > 1500) template spectrum for each of our target 
stars by combining all the individual dithered spectra using an Interactive Data Language (IDL) 
routine which performs neighborhood averaging on unevenly sampled data. This template spectrum 
is used to estimate bisector errors (§2.2). 
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2.1. Line Bisector Analysis 

We construct the bisector for the Fe I 6252.6 A line following the method described by Gray 
(1988). Specifically, we locate the pixel closest to the line core and then step up the blue side of 
the profile pixel by pixel. From each blue pixel, we extend a horizontal (i.e., same level of residual 
intensity) line segment across the profile. The wavelength value for the red end of the segment is 
determined though linear interpolation between the neighboring pixels. The bisector connects the 
midpoints of these line segments. We adopt linear interpolation because error propagation is more 
straightforward, bisector points are more independent of one another, and outliers do not cause 
ringing. There is likely some curvature across a pixel, especially in the line core, but we fitted the 
line core with a parabola and curvature is minimal in the steep part of the line profiles, where we use 
linear interpolation. In addition, the points on one side of the profile are used without interpolation 
which further mitigates the effect of curvature. Typically, we used seven bisector points for Fe I 
6252.6 A in G-type stars and twelve bisector points for the more rotationally broadened F-type 
stars. A typical line bisector is illustrated in Fig. 2. 

Once we have derived the bisector, we measure velocity offsets vi, V 2 , and V 3 , corresponding 
to residual intensities yi, 1/2, and 1/3 (see Fig. 2). Velocity offsets are measured with respect to the 
center of a parabolic fit to the line core. We determine V{ for a given y, by linear interpolation of 
adjacent bisector points. Following Gray (1998) and Gray & Hatzes (1997), we chose ( 1 / 1 , t/ 2 , 1 / 3 ) = 
(0.83,0.7,0.5) for the Fe I line, which for 51 Peg has a residual intensity of ~ 0.3 in the core. For 
other stars, the measurement depths, y, were scaled according to observed line depth. 

Line asymmetries are typically characterized by velocity span (or amplitude) or bisector cur- 
vature. For example, the former approach is used by Livingston et al. (1999) to measure line profile 
asymmetry in solar spectra. Livingston et al. (1999) define the bisector amplitude as aj, — A c 
where A c is the observed central wavelength of the absorption feature and x mtn is the most extreme 
blue point in the bisector. This approach is especially appropriate for solar observations where the 
combination of very high spectral resolutions (R > 3 X 10 5 ) and S/N (> 5000) yields well-defined 
bisectors with a large number of independent points. Consequently, derivatives of the bisector itself 
can be determined and x min can be deduced through minimization. This method is usually not 
feasible for stellar observations which are typically characterized by lower spectral resolution and 
S/N ratio. In the context of this investigation, we define the velocity span as 


v s = v 3 -v 1 . (2) 

The other principal measure of line asymmetry is bisector curvature, which has been described 
as the difference between bisector spans computed for the top and bottom halves of a line (Gray 
& Hatzes 1997; Gray 1998). In our terminology, this is equivalent to 


Cb = («3 ” u 2 ) - (t>2 - ui)- 


( 3 ) 
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Since the curvature, c&, is a second derivative (whereas span is a first derivative), the signal- 
to-noise requirements are more demanding in order to obtain significant measures of variability in 
the bisector curvature. We therefore adopted an alternative characterization which is a variant of 
the span method designed to increase sensitivity. We define this velocity displacement as 


tq + v 2 + v 3 . 

Vb = g A c , 

where Vb is the average displacement of the evaluation points from the line core. 


( 4 ) 


A comparison of the formulations in equations (2)-(4) suggests that for a “C”-shaped bisector, 
Vb and v s will be negative while the curvature, q,, will be positive. The principal advantage of Vb 
as defined by equation (4) is that by utilizing the average of all three evaluation points, the impact 
of the uncertainties in individual evaluation points is mitigated. We also note that our method 
of bisector measurement relies on an accurate determination of the location of the line core (A c ), 
which we ascertain by fitting a parabola to the line core. 


In equation (4), averaging tq, v 2 , and v 3 reduces uncertainties relative to the separate use 
of these points in a first or second derivative. Uncertainties in the velocity shift of the line core 
are reduced by fitting a parabola to a few points in the line core. In essence, we are utilizing 
two averaged points for the determination of the bisector amplitude. In contrast, the curvature 
formulation in equation (3) uses three independent points with one point (v 2 ) used twice. We find 
from our data that this latter approach significantly increases measurement uncertainties (see § 
2 . 2 ). 

Finally, we do not include any additional bisector points in the construction of Vb since shallower 
points are intrinsically noisy and only further increase the uncertainty in Deeper points impinge 
on the line core, which would compromise the independence of our two bisector span reference 
points. 


2.2. Error analysis 

To estimate noise in the continuum, we divide each spectrum by a summed template spectrum, 
which is essentially noiseless. We then calculate the standard deviation of this ratio in continuum 
windows where the residual intensity of the template is between 0.95 and 1.02. We assign un- 
certainties (a 5 ) to every point in the spectrum by multiplying the continuum uncertainty by the 
square root of the residual intensity. In other words, we assume that uncertainties are controlled 
by Poisson statistics in the stellar signal, which is approximately true for these bright stars. 

Uncertainties in residual intensity first propagate into our bisector determination when in- 
terpolating the red wing of the profile as illustrated in Figure 2. Uncertainty in the wavelength 
corresponding to a particular residual intensity is given by 
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< 7 w = a s jm , (5) 

where m is the slope of the spectrum (Gray 1988), which is measured using the essentially noiseless 
template. We then propagate errors through successive steps in the bisector construction and 
interpolation onto reference depths (y) to obtain Oi, a 2 , and a 3 , which are the uncertainties in 
uj, v 2 , and v 3 respectively. Uncertainty in the wavelength of the line core, a\ c , is derived from 
estimated uncertainties in the residual intensity points being fitted. The final error in v b is 


_2 _ °\ + °\ + a 3 , 2 

*<■ ~ 9 + ^ 


( 6 ) 


The observed scatter in bisector velocity displacement is consistent with the formal uncertainties 
calculated in this manner, as discussed below. 


We also examined the uncertainties associated with other methods of characterizing bisectors, 
namely span and curvature. In particular, the uncertainty in span from equation (2) is: 


= ol + ul, 

while the uncertainty in curvature from equation (3) is 


( 7 ) 


a] = a\ + Aa\ + a^. (8) 

In Table 2, we compare various characterizations of line bisectors for 51 Peg, namely velocity 
displacement (uj,), span (u s ), and curvature (c;,). For each technique, we give the weighted mean 
of all measured values for 51 Peg (column 2), the standard deviation about the mean (column 3), 
and the average formal uncertainty (column 4). The last column contains the observed standard 
deviation divided by the average formal error. Values near unity indicate that our estimates of 
formal uncertainty provide a realistic indicator of actual errors. 

Table 2 shows that for our data, velocity displacement as defined in equation (4) yields smaller 
formal errors than span or curvature. Moreover, the mean velocity displacement is as large as 
other characterizations, making v b the most sensitive diagnostic of variability. Since our purpose 
here is to search for periodicity in bisector behavior, we adopt the formulation that provides the 
strongest constraints on bisector amplitude variations, namely velocity displacement. For illustra- 
tive purposes only, we present in Figure 3 bisector curvature and span versus orbital phase for 51 
Peg. 
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3. Results and Discussion 
3.1. Significance of Variability 

In Table 3 we summarize our measurements and formal errors. Column (4) gives the mean 
of vt, for individual observations, weighted by their corresponding uncertainties. Column (6) gives 
the formal uncertainty in the weighted mean, calculated from uncertainties in Vb for each 
observation. See equation 4.31 of Bevington & Robinson (1992) for details. Column (5) gives the 
mean value of the formal errors, (<r). Measurements for each night are depicted graphically as time 
series in Fig. 4. For convenient reference, (t>b), 07, and (a) are given in the figure. 

A key goal of this investigation is to test whether line asymmetries vary periodically on time 
scales consistent with precisely measured radial velocity variations. For each star, we plot in Fig. 5 
bisector velocities as a function of Julian date and orbital phase, as computed from planetary orbital 
periods and reference times (T 0 ) in Table 1. At the reference time, the radial velocity transitions 
from negative to positive. Inspection of Fig. 5 reveals no obvious periodic variability in bisector 
velocity displacement, though excursions for individual observations may be marginally significant 
in some cases. 

In Table 4 we present the ratio of the observed scatter, a v , to the mean of the formal errors, 
{a). This quantity is analogous to the last entry in the first row of Table 1, but for all stars. For 
ratios greater than unity, bisector variations may be significant, even if not periodic. A comparison 
of the mean value of this ratio for stars with planets compared to those without planets (10 Tau, 7 
Ser, and 16 Cyg A) yields values of 1.1 and 0.9, respectively. We conclude that there is no significant 
statistical difference in the behavior of line asymmetry in stars with and without detected planets. 


3.2. Periodogram Analysis 

An important question is whether periodic variations in bisector amplitude translate into 
periodic, apparent Doppler shifts in line profiles observed at moderate resolution, thus mimicking 
the radial velocity signature of a very low mass (planetary) companion. While it is plausible that 
variations in line asymmetry may lead to subtle, apparent Doppler shifts, quantifying this assertion 
is difficult. In the absence of suitable hydrodynamic calculations or even empirical data relating 
changes in line bisector amplitude to apparent Doppler shifts, we can only crudely estimate the 
magnitude of the effect. Certainly, we can expect that any apparent change in stellar radial velocity, 
Ai) r , would be such that Av r < Ad*,. Inspection of Fig. 1 suggests that Av r < Au(,/2 is a more 
realistic estimate of the magnitude of the change. Thus, changes in bisector amplitude of at least 
100 m s -1 would be required to create intrinsic variations comparable to the 56 m s -1 radial velocity 
variations seen in 51 Peg. 

We performed a periodogram analyses of Vb time series for each star (Fig. 4), using an algorithm 
from Press &: Rybicki (1989). Upon computing false alarm probability levels and examining the 
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periodograms, we did not find any significant peaks at the planetary orbital period or any other 
period. This test is less sensitive than specifically searching for periodic variability at the period 
determined by precision radial velocity studies, but provides new constraints on the existence of 
periodic variations on other time scales. 


3.3. Sensitivity Estimates 

We also ran simulations incorporating our measurement uncertainties and actual observing 
windows to determine our sensitivity to signals varying at the planetary orbital period. The results 
indicate that except for v And, we could have detected the presence of a continuous periodic signal 
at the planetary orbital period with an amplitude equal to the measured amplitude of radial velocity 
variations. For v And, errors in measured bisector amplitude are too large to permit detection of 
a periodic signal with this amplitude. The high vsini and the relatively low amplitude of radial 
velocity variation, along with our errors, would mask any regular variations in vi, that might be 
present. In all other program stars with planets, the periodogram readily revealed the simulated 
periodic signal with false alarm probabilities that ranged from ~ several X 10 -2 to ~ 10 -11 . We 
conclude that for these stars, changes in line asymmetry are not the major cause of observed radial 
velocity variations. 

For illustrative purposes, we present in Fig. 6 the periodogram for the 51 Peg time series of 
observations along with a simulation that includes a periodic signal at the orbital period and the 
velocity amplitude of the reported planet. In the specific case of 51 Peg, we performed simulations 
for a range of amplitudes extending from 5 m s -1 to 50 m s” 1 . The results reveal that variations 
at the orbital period with an amplitude of only 25 m s" 1 are readily detectable within our errors 
at a false alarm probability ~ 10“ 3 . But at amplitudes as small as 10 m s -1 , any such regular 
variations that may be present would not be discernible in our data. 


3.4. Trend with Spectral Type 

Finally, we display in Fig. 7 mean bisector velocity displacement (t>&) as a function of B — V 
color. Inspection of Fig. 7 suggests that the { v 5) is correlated with stellar spectral type, a result 
that is consistent with Gray (1988). In particular, Gray (1988) notes that the redward turn of the 
bisector in the line core (the bottom part of the “C”) is absent in early-type stars, only manifesting 
itself near spectral type GO. We would expect stars showing only the top of the “C” to have 
positive and this is indeed what we observe for most of our F stars. For cooler stars, the 
velocity displacement becomes more negative. Note that 55 Cnc (G8V) is the only star in our 
sample with molecular CN lines visible in our summed templates. The blue wing of the Fe I 6252.6 
A line is blended with a line of CN, which may affect our bisector characterization for 55 Cnc. 

The significant departure of r Boo from the trend shown in Fig. 7 may be related to its 
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relatively rapid rotation period of 3.3 i 0.5 days (Baliunas et al. 1997). As discussed by Gray (1988), 
and confirmed in the case of the Sun by Smith, Huang, & Livingston (1987), line asymmetries as 
characterized by the amplitudes of line bisectors can be enhanced by rotational redistribution of the 
Doppler shifts arising from the granular motions. In essence, rotation adds a velocity component 
to the line-of-sight motion of granules with the effect most pronounced at the equatorial limbs. 
Certainly, the results in Fig. 7 for r Boo are qualitatively consistent with the nature of this 
so-called “rotation effect.” 

Alternatively, the coincidence of the stellar rotation period with the 3.3 day orbital period of 
the planetary companion raises the speculation that the rotational synchronization via tidal locking 
of the planet with the star could lead to the modification of the granulation pattern in the pho- 
tosphere. Cuntz, Saar, & Musielak (2000) conclude from theoretical and empirical considerations 
that extrasolar giant planet systems with small separations such as r Boo should be characterized 
by both strong magnetic and tidal interactions. In fact, r Boo exhibited the strongest magnetic 
interaction among the stars considered in their theoretical investigation. In view of the well-known 
interaction between magnetic fields and convection in the Sun and late-type stars, this speculation 
merits investigation. 


4. Conclusions 

We have utilized synoptic spectroscopic observations of stars with reported planetary compan- 
ions to investigate whether line profile variability arising from stellar atmospheric motions could 
mimic the observed radial velocity variations. In summary, we did not find any periodic variability 
in the velocity amplitude of the bisector of the photospheric Fe I line at 6252.6 A which was the 
focus of our investigation. Furthermore, the variations in bisector amplitudes that were detected 
are consistent with uncertainties expected from the propagation of errors in the original spectra. 
While isolated excursions in the time series of bisector amplitudes are observed that may arise from 
real variability, the significance is not high enough to exclude other sources of measurement error. 

The behavior of our measured mean bisector velocity displacement as a function of spectral 
type is consistent with previous work. While line profile variations are certainly a source of intrinsic 
stellar noise in Doppler spectroscopic searches at the sub-km s" 1 level, our results indicate that it is 
the presence of substellar companions which best explains the periodic variations in radial velocity 
that have been reported for our program stars. 

We gratefully acknowledge support by a grant from NASA to the NOAO/NSO under the 
auspices of the Origins of Solar Systems Program which made this investigation possible. 
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Fig. 1. — a-d. Template spectra of all the stars in our data set. Spectral resolution (A/JA) is ~ 1.2 x 10 5 
with S/N> 1500. The F-type stars (see Table 1) exhibit shallower, rotationally broadened absorption lines. 
The spectrum of 55 Cnc contains many more lines than are seen in our other target stars, most likely due 
to an enhanced CN content in the photosphere. 
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Fig. 2. — Schematic diagram of a spectral line bisector, illustrating the “C”-shape and evaluation 
points ui,u 2 , and v 3 (plotted as x). The bisector points (small, filled circles) are the midpoints 
of the horizontal, dashed-dotted line segments, plotted on the wavelength scale of the spectrum. 
Filled squares are data points in the observed spectrum, while open squares indicate where linear 
interpolation has been used to determine the red endpoints of the bisector segments. The solid line 
with large, open circles is the entire bisector, plotted on the velocity (v) scale, which is expanded 
by a factor of 40 relative to the wavelength scale. Zero velocity is defined as A c , determined from 
a parabola fit to the line core (dashed curve). 
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Fig. 4. — Time series of bisector amplitudes for the program stars. The dotted line on each plot 
denotes the weighted mean value of v i calculated using the formal errors. Also given are a v , 
the variance of the distribution of Vf, and (a), the mean value of the formal errors. 
















Fig. 5. — Plots of Vb versus phase <j) for those stars exhibiting periodic Doppler velocity variations 
due to the presence of a substellar companion. Orbital phase of the companion was computed from 
the parameters listed in Table 1. 
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Fig. 6. — Periodograms computed for (a) the observed time series for 51 Peg and (b) a simula- 
tion assuming a periodic signal at the planetary orbital period. The simulation incorporates the 
observing window and errors that are present in (a). No significant peaks are seen in the obser- 
vational data. A significant feature is evident in the simulation at a frequency of 0.2373 day -1 , 
corresponding to a period of 4.214 days, near the planetary orbital period of 4.231 days. 
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B - V (mag) 

Fig. 7. — Mean bisector velocity displacement plotted against spectral type of each star. The 
transition from negative to positive values of (i^) towards earlier-type stars could be interpreted as 
movement up the “C” (see text). 
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Table 1. Target Stars 


Period a K 

HR Name Type (days) T 0 a?6 (m s -1 ) 


458 

v And c 

F8 V 

1101 

10 Tau 

F9 IV-V 

3522 

55 Cnc 

G8 V 

4277 

47 UMa 

Gl V 

5072 

70 Vir 

G4 V 

5185 

r Boo 

F6 IV 

5933 

7 Ser 

F6 IV-V 

5968 

p CrB 

GO Va 

7503 

16 Cyg A 

G1.5 Vb 

8729 

51 Peg 

G2.5 IVa 


4.6163 

50542.736 

60.1 

14.6597 

50484.188 

72.8 

1026.58 

50405.303 

47.1 

116.666 

50539.287 

317.0 

3.31253 

50537.569 

468.2 

39.8412 

50975.357 

60.7 

4.23110 

50509.695 

56.1 


Calculated from data courtesy G. Marcy. 
b JD— 2400000 

c Multiple-planet system — Doppler parameters are given for 
the innermost companion only (v And b). 


Table 2. Bisector Amplitudes for 51 Peg 


Measure a 

Mean 
(m s” 1 ) 

Scatter 
(m s" 1 ) 

Error 
(m s" 1 ) 

Ratio 

Vb 

-15 

42 

29 

1.45 

V s 

2 

83 

61 

1.36 

Cb 

14 

53 

88 

0.60 


a Bisector amplitudes: v b is defined by equation 
(4); v s is given by equation (2); q, is calculated ac- 
cording to equation (3). 
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Table 3. Summary of Results 


HR Name Type (wj,) (<r) a ^ N 0 (, s 

(m s -1 ) (m s _1 ) (m s _1 ) 


458 

v And 

F8 V 

1101 

10 Tau 

F9 IV-V 

3522 

55 Cnc 

G8 V 

4277 

47 UMa 

Gl V 

5072 

70 Vir 

G4 V 

5185 

r Boo 

F6 IV 

5933 

7 Ser 

F6 IV-V 

5968 

p CrB 

GO Va 

7503 

16 Cyg A 

G1.5 Vb 

8729 

51 Peg 

G2.5 IVa 


12 

81 

16 

64 

31 

36 

6 

80 

-41 

36 

4 

225 

-3 

43 

6 

121 

-17 

32 

4 

135 

-45 

200 

22 

118 

46 

129 

15 

120 

5 

33 

4 

105 

-15 

42 

5 

82 

-15 

29 

4 

115 


Table 4. Ratios of Scatter to Mean Uncertainties 


HR Name Type ct„/(<t) 


458 

v And 

1101 

10 Tau 

3522 

55 Cnc 

4277 

47 UMa 

5072 

70 Vir 

5185 

T Boo 

5933 

7 Ser 

5968 

p CrB 

7503 

16 Cyg A 

8729 

51 Peg 


F8 V 

0.7 

F9 IV-V 

1.1 

G8 V 

1.5 

Gl V 

1.1 

G4 V 

1.2 

F6 IV 

0.7 

F6 IV-V 

0.6 

GO Va 

1.4 

G1.5 Vb 

0.9 

G2.5 IVa 

1.4 




